Quantum adiabatic polarons by translationally invariant perturbation theory 



0. S. Barisic 1 and S. Barisic 2 
1 Institute of Physics, Bijenicka c. 46, HR-10000 Zagreb, Croatia 
2 Department of Physics, Faculty of Science, University of Zagreb, Bijenicka c. 32, HR-10000 Zagreb, Croatia 

The translationally invariant diagrammatic quantum perturbation theory (TPT) is applied to the 
polaron problem on the ID lattice, modeled through the Holstein Hamiltonian with the phonon fre- 
quency ljo, the electron hopping t and the electron-phonon coupling constant g. The self-energy dia- 
grams of the fourth-order in g are calculated exactly for an intermittently added electron, in addition 
to the previously known second-order term. The corresponding quadratic and quartic corrections 
to the polaron ground state energy become comparable at t/uio > 1 for g/uio ~ (t/uo) 1 when the 
electron self-trapping and translation become adiabatic. The corresponding non adiabatic/adiabatic 
crossover occurs while the polaron width is large, i.e. the lattice coarsening negligible. This result is 
extended to the range (t/ujo) 1 ^ 2 > g/wa > (t/cuo) 1 ^ 4 > 1 by considering the scaling properties of the 
high-order self-energy diagrams. It is shown that the polaron ground state energy, its width and the 
effective mass agree with the results found traditionally from the broken symmetry side, kinematic 
corrections included. The Landau self trapping of the electron in the classic self-consistent, localized 
displacement potential, the restoration of the translational symmetry by the classic translational 
Goldstone mode and the quantization of the polaronic translational coordinate are thus all encom- 
passed by a quantum theory which is translationally invariant from the outset. This represents 
the first example, open to various generalizations, of the capability of TPT to hold through the 
adiabatic symmetry breaking crossover. Plausible arguments are also given that TPT can describe 
the g/wo > (t/cuo) 1 ^ 2 regime of the small polaron with adiabatic or non-adiabatic translation, i.e., 
that TPT can cover the whole g/u^o, t/tvo parameter space of the Holstein Hamiltonian. 

PACS numbers: 71.38.-k, 63.20.Kr 



I. INTRODUCTION 

The polaron is one of the earliest examples of a topolog- 
ical particle associated with symmetry breaking. Similar 
examples in the condensed physics are phasons soli- 
tons 21, magnetic polarons [Sh Kosterliz-Thouless vor- 
tices 1|) Zhang- Rice singlets @, and manyothers. Po- 
laron was initially introduced by Landau jfj as a state 
of an electron coupled to the classical deformation field, 
which, by adiabatic self-localization of the electron (elec- 
tron self-trapping) , breaks the original translational sym- 
metry of the Hamiltonian. The corresponding symmetry 
restoring Goldstone mode (classical or quantum) is then 
the translation of the polaron 

This description of the (adiabatic) polaron, from the 
broken symmetry side, suggests that the translationally 
invariant quantum perturbation theory (TPT) in terms 
of the electron-phonon (e-p) coupling constant g, can- 
not reach such a (symmetry restored) state. Indeed, the 
symmetry breaking (at T = 0) is, as a rule, associated 
with the singularity (quantum critical point, QCP) in the 
ground state energy as a function of g. This is usually 
taken to restrict the use of TPT to the high symme- 
try (translationally invariant) phase on one side of the 
QCP. On the other hand, the study of some continuous 
e-p Hamiltonians (including the 3D Frohlich Hamiltonian 
which exhibits a polaron as a classic solution) , concluded 
that the correspondi ng q uantum ground state energy is a 
smooth function of g This shows that, depending on 
the dimension of the system and the range of the forces, 
the quantum fluctuations of the deformation field can re- 



move the QCP and replace it by a smooth crossover. 
However, this does not gua rantee that TPT can hold 
through such a crossover (lulL2ll3l |. because the absence 
of a singularity for g's, which make the Hamiltonian her- 
mitian, can say nothing about the radius of convergence 
and the behavior of the perturbation series in the com- 
plex g plane. The question is not only whether TPT can 
hold through the symmetry breaking crossover in prin- 
ciple 0| but also how and for which physical regimes 
151 this can be done in practice. Some elements of the 
answer to this question are given here. 

In the adiabatic limit, which is of the main inter- 
est here, the electron is "self-trapped" (adiabatic self- 
trapping), staying always in the same localized state, 
which moves with the polaron distortion of the lattice. 
In the continuous approximation, the adiabatic polaron 
is free to move along the lattice. However, the lattice 
coarsening introduces the Peierls-Nabarro (PN) or Umk- 
lapp potential 0, 0] into the motion of the adiabatic 
polaron. If treated classically, even a tiny periodic poten- 
tial pins the polaron of mass M p to the lattice. There are 
thus two possible, separate, mechanisms which can break 
the translational symmetry of the lattice on the adiabatic 
level. The first is the adiabatic self-trapping of the elec- 
tron. The second is the pinning of the adiabatic polaron 
to the PN potential of the lattice. However, although 
there are two symmetry breaking mechanisms, there is 
only one symmetry to break, namely the translational 
symmetry of the lattice. It is therefore important to 
distinguish between the electron adiabatic self-trapping 
and the polaron pinning. As the electron adiabatic self- 
trapping is prerequisite for the polaron pinning, it is the 
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fundamental mechanism of symmetry breaking. For this 
reason, the problem of the applicability of the TPT is fo- 
cused here on the electron adiabatic self-trapping on the 
lattice. 

In fact, the inter-relation of the electron self-trapping 
and the polaron pinning can be fully understood in the 
adiabatic regime. As illustrated schematically in Fig. 
this regime can be reached from the weak and from the 
strong coupling side. For small g, the quantum electron- 
phonon system is described reasonably well by the low- 
est order perturbation theory in the e-p coupling, which 
is apparently translationally invariant and nonadiabatic 
[23| • F° r the discrete lattice the polaron motion is nona- 
diabatic also for extremely large g. The electron must 
use phonons nonadiabatically to leave its (too) strongly 
pinned adiabatic phonon correlation cloud in order to 
gain the derealization energy. Obviously, this quantum 
state is also translationally invariant. The adiabatic po- 
laron regime, if it exists, is thus separated from the two 
nonadiabatic, translationally invariant regimes either by 
a pair of QCP's or by a pair of corresponding crossovers. 
Which is the case can be found from either (small or large 
g) side. Again however, in the case when QCP's are re- 
placed by crossovers, it remains to be proven that TPT 
can hold through them. 

To make the picture complete, it should be further re- 
alized that on coming from the low g side, the QCP or 
the crossover corresponding to the electron adiabatic self- 
trapping can, generally speaking, occur either towards 
the large adiabatic polaron quantum state, when the PN 
potential although present is negligible, or towards the 
small adiabatic polaron quantum state, when this poten- 
tial plays an essential role. In the former case there is an 
additional crossover, the one from large to the small adia- 
batic polaron, as indicated by the dashed curve in Fig. 
Importantly, this behavior can never correspond to the 
breakdown of the translational symmetry of the lattice 
(i.e., to QCP), because, as already mentioned, the latter 
is consumed by the electron adiabatic self-trapping QCP 
or by the corresponding crossover. It appears therefore 
that the electron adiabatic self-trapping on the lattice 
provides a critical test for the applicability of TPT. In 
this context it is natural to study the low g side of the 
problem, because TPT is the expansion in terms of g. 
If TPT is able to hold through this crossover or QCP it 
is likely that it applies, at least in principle, to all val- 
ues of g. But then, so does also the theory which starts 
from the broken symmetry side, and it is only a matter 
of convenience which approach is to be used when. 

In the present paper, TPT is tested on the discrete 
ID Holstein model. This model provides a simple ex- 
ample of the general situation described above. For 
sufficiently fast electrons it exhibits a pair of nonadi- 
abatic/adiabatic crossovers, which replace QCP's [13]. 
The small g crossover corresponds to the electron self- 
trapping in the state of the large adiabatic Holstein po- 
laron 0, |2i| ■ At sufficiently large g it is supplemented 
by the large / small adiabatic polaron crossover in the PN 
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Figure 1: Schematic phase diagram for the polaron problem 
in terms of the e-p coupling constant g and the electron mass 
m e i (arbitrary units). Full lines denote QCP's or crossovers, 
which could be responsible for the breaking of the transla- 
tional symmetry of the lattice. The broken line denotes the 
crossover that is associated to the pinning of the polaron to 
the lattice. 



potential, which is thus well separated from the electron 
self-trapping crossover. It will be argued here that TPT 
can reach beyond the low g electron self-trapping show- 
ing explicitly that this QCP is replaced by the crossover 
and that the lattice coarsening (PN, Umklapp) effects are 
negligible in this case. Finally, some plausible arguments 
are given which indicate that TPT is valid for all val- 
ues of g, although the full proof of this assertion requires 
additional considerations. 



II. GENERAL 

The Holstein Hamiltonian on the discrete ID lattice 
with L sites is given by (H= 1) 



L/2-1 

-L/2 



-t 4(c r+ i + Cr-i) + b\b r - g clc r (bl + b r )} 



(1) 

where the fermion and boson operators c r and b r on the 
site r are defined in the usual way. Equation is meant 
to describe N electrons subjected to hopping t along the 
chain and to the local interaction g with the displace- 
ments of the lattice u r = Xo(bj. + b r ). At g — the latter 
behave as harmonic oscillators with frequ ency uo and the 
zero-point displacement xq = y/l/2MiJo, where M is the 
ionic mass. The properties of the Holstein Hamiltonian 
for a given L and N are thus described in the simple 2D 
parameter space, e.g. in terms of g/ui and a = uj /2t. 
The only adimensional quantity independent of M that 
can be constructed from g/ojQ and a is A = a g 2 /ujQ, 




Figure 2: The exact electron self-energy E' '. The double line 
is the exact G, the shaded triangle is the exact vertex and the 
wavy line is the free phonon propagator Do- 



which has to appear naturally (instead of g 2 /ojq) in the 
adiabatic regime, of interest here. 

The structure of TPT for the Hamiltonian Q and sim- 
ilar models was examined in Refs. [22L I23L l23 | . The main 
conclusion in Ref. [24| was that the polaron properties 
can be conveniently extracted from the correlation func- 
tions which form the basis of the diagrammatic version 
[2H of TPT. Two types of correlation functions are to be 
distinguished in this respect, namely those which by con- 
struction conserve the number N of electrons, and those 
which do not. The displacement-displacement correla- 
tion function D belongs to the first class, and the electron 
propagator G to the second, because it changes N by ±1. 
In the polaron case D = D^ 1 ' is to be taken at N = 1 and, 
consistently, G = G^ at N = 0. G' ' describes the prop- 
agation of an electron intermittently added to the system 
of bosons. Actually, on using the Lehmann representa- 
tion of G^ it turned out j2J| that the most important 
properties of the polaronic correlations can be calculated 
from the corresponding irreducible electron self-energy 
alone. E^ ) gives the position of the polaron bands, 
in particular the value of the ground state energy, the ef- 
fective mass of the polaron, and, in the continuous limit, 
the polaron width, as will be further discussed below. 

The reduction of the problem to the calculation of E^ ) 
represents an important simplification, because at N = 
all fermionic (Pauli) correlations are eliminated from the 
outset. This contrasts with the calculation [24[ of the 
ground state energy from its direct diagrammatic expan- 
sion at N = 1, where the electron exchange effects dis- 
appear from the result only after tedious cancellations. 

This important simplification appears formally in the 
exact expression for given in Fig. [2] in the usual 

[2ll ] diagrammatic language. In contrast to the general 
N case, the wavy line in Fig. [^represents the free N = 
displacement-displacement correlation function = 
A), 
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Indeed, at N = 0, the phonon renormalization, which 
starts necessarily with the creation of the electron-hole 
pair, is impossible: there is no electron in the system, 
additional to that created intermittently. In other words, 
the electron described by G^ can only advance in time 
in Fig. |2 until it is annihilated, i.e., G^ can have poles 






Figure 3: Second- and fourth-order diagrams for E^. 

only in the lower w-half-plane. This holds in particular 
for the free electron propagator G^ , 



G { °\k,co) = 



1 



w - & + if) 



(2) 



For further convenience the zero of the free electron en- 
ergy £k in Eq. <j2j is taken, unlike in Eq. QJ, at k = 0, 



£ fe = 2t(l-cos(fc)) . 



(3) 



Figure also contains the triangular vertex correc- 
tion of Fig. originally discussed by Migdal in 
the limit of a large number of electrons N. At large 
N, especially when the soft-phonon renormalizations of 
the phonon propagator are important, this vertex cor- 
rection can sometimes be neglected. However, at N = 0, 
when D ( °) = D Q in Fig. the vertex correction can be 
as important as the corresponding self-energy correction, 
as will be further seen below. The fact that the "soft- 
phonon" corrections do not appear at N — in Fig. |2] 
does not mean that the phonon propagator at N = 1 
is not developing a soft-phonon branch, as a signature of 
the polaronic correlations. It does, but this branch does 
not enter the calculations of E' ' in the small- N hierar- 
chy of the correlation functions. 

In order to emphasize that at = the electron can 
only advance in time, the irreducible E^ ' of Fig. |2]can 
be represented by the expansion, order by order in g. In 
Fig. the arrow on the free electron line G^ of Eq. 
points explicitly forwards in time. As usually, the 
diagrams in Fig. represent the perturbation series in g 
either in the direct or in the reciprocal space. Concerning 
the propagation in space it should be noted that E2 and 
the non-crossing contribution T,^ c are local because so 
is D of Eq. J31 , in contrast to the crossing diagram Ef 
which contains also the phonon-assisted intersite prop- 
agation. In terms of Fig. [21 E4 7 represents the leading 
(Migdal) vertex correction to E^ ' . It will be shown below 
that in the large adiabatic polaron regime the non-local 
contributions of the crossing diagrams to the low-energy 
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properties are equally important for the polaron forma- 
tion as the local corrections of the NC diagrams. 

The low-energy properties are associated with the po- 
laron spectrum Assuming that the perturbation series 
for S(°) shown in Fig. El is meaningful, is the solution 
of the equation 



with 



9^ L/2 ~ 1 

'<•> = X E 



/ e + 2 sin 2 (ttto/L) 



(8) 



I* =&+iteE (0) (&,&), (4) 

where ReT,^ represents the real part of E^ ). It is im- 
plicit in the concept of the quantum polaron, as of the 
coherent quantum state, that the imaginary part of E*- -* 
is inflnitesimally small around the polaron pole This 
is shown later to be consistent with the present analysis 
and, anticipating this result, no distinction is made fur- 
ther between ReT,^ and E^ itself. The energy of the 
k = state in the lowest polaron band is obviously re- 
lated to the ground state energy. Noting that the ground 
state energies Eq (0) and Eq ^ (0) of the noninteracting 
(g = 0) N = 1 system and the N = system, respec- 
tively, are taken to coincide by choosing £/_o = in Eq. 
<EJ , the ground state energy gain A of the N = 1 system 
with finite g is 

A = (0) - EP (g) = -E(°) (k = 0, - A) , (5) 

according to the Lehmann representation of . It is 
worthy of note that the higher energy solutions of Eq. 

give the values of the polaron k = energies of the 
excited polaron bands. It is reiterated that the upper 
index in Eq. (JjJ and elsewhere denotes the number of 
electrons involved. 

The behavior of £fc in the vicinity of k = defines 
the polaron mass M p which is given by the derivates of 
E(°)(£ fc ,u;) at k = and w = - A, 



where 



M" 



2t 



1 + d£(°> /dj k 
l-dE(°)/dw 



(6) 



It is worth noting that the denominator in Eq. © taken 
at arbitrary k gives the corresponding electron spectral 
density. 

Finally, it will be shown latter that the polaron width 
can be related to A of Eq. (0 , in particular for the large 
adiabatic polaron. This completes the list of the main 
polaron properties (energy, mass, width) which can be 
derived from E^ ) (k, uf). The next step is the examination 
of £<°>(fc,w) itself. 



III. LOW-ORDER DIAGRAMS 

The behavior of the leading (skeleton diagram) dia- 
gram E2 (dropping henceforth the (0) superscript) of Fig. 
El is quite instructive J24[. E2 is given as 
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and the integer m is related to the phonon wavevector q 
by the Born- van Karman boundary conditions, 



q = 2Trm/L . 



(10) 



The sum © can be expressed in the closed, analytical 
form for arbitrary e > 0. The regime of interest here is 
L~ 2 < e <C 1. For e > L~ 2 finite size effects can be 
neglected, i.e., the sum in Eq. (jHJ can be turned into 
the integral over q. If e <C 1, q c = \[2e plays the role 
of the infrared cutoff in this integration. In other words, 
2 sin 2 (q/2) = 1 — cos(g) in Eq. JHJ can be replaced then 
by q 2 /2 and integrated from q c to infinity, yielding 



(11) 



The same result is obtained by taking the continuous 
limit in the Hamiltonian Q from the outset, except that 
the present procedure shows that the corrections due to 
the lattice coarsening (to the Umklapp processes) are of 
the order of e <?C 1 itself. 

In the weak-coupling limit, when the ground state en- 
ergy shift A of Eq. © is small A < wo, w = -A can be 
neglected in e of Eq. ^ appearing in combined equations 
© and l(TT|l . Taking thus e ~ a <C 1 leads to 



wo 



2tAy^ 



(12) 



This apparently is a nonadiabatic result, because A is 
dependent on the ionic mass M through a. In the weak- 
coupling limit L~ 2 < e « a -c 1 thus ensures vali dity 
of the continuous [24[ rather than of the adiabatic (22 | 
approximation. Alternatively, the inequality L~ 2 < e <C 
1 may be thought |24[ of as fixing the order of the limits, 
first L — > 00, and then t/wo — > 00, in the search for a 
nontrivial adiabatic regime. 

In this context, the nonadiabatic nature of Eq. I|12|) 
for a < 1 suggests the calculation of the higher order 
diagrams E4 and £4 of Fig. 0in the first place. The 
two internal frequency integrations in those diagrams can 
be carried out easily to give 



5 



S^ c = JL_/'( e )/( e + a) 



(2t) 3 



(13) 



sy = 



(20 3 ^ [e + 2 sin 2 (g/2)] [e + 2 sin 2 (q'/2)} 

(14) 



[£ + a + 2sin^((/c + g-g')/2)] 

where /(e) in Eq. (TT3")) is given by Eq. {jJJ, /'(e) is its 
derivative with respect to e, while q, q' and k in Eq. <14ll 
are given by Eq. tflTTfl . Translational invariance was used 
in Eq. ifH)) to associate the dependence of £4 on the 
external momentum k with the convolution of q and q', 
i.e., with the crossing of two phonon lines in Fig. In 
contrast, as already mentioned, £4 (like £2) is local, 
i.e., independent of k. 

While /(e) and thus I'(e) of Eq. l|T3"j) are known in 
the closed form for e > arbitrary, £4" calculated in 
Appendix is exhibited here only in the interesting limit 
L~ 2 < e + a <C 1, when the q, q' summations in Eq. 
Ill4t can be turned into the infrared singular integrations. 
Together with Eqs. ltl3|l and ifTTjl this gives 



(2i)3 ^[eT, 



2^e + Ve + a 



(15) 



(2t)3 e^FT^ k 2 + 2(2^ + VF+otj 2 

both results holding with the accuracy l/(e + a). 

Equations ltl5)l and iflfijl can be used first to define the 
range of validity of the low-order perturbation theory by 
comparing £2 and £4 at k = and lo = — A of Eq. (JJJ. 
The two become comparable for A comparable to ljq, i.e., 
for 



polaron mass, given by Eqs. © and respectively. 
This regime is thus beyond the reach of the DMFT [22]. 

As has been already emphasized before, the condition 
(tl7l when introduced in Eq. ltl2)l leads to 



2iA 2 



(18) 



for a <C 1, where, once again, A = g 2 /2toJo is indepen- 
dent of ionic mass. The estimate ijlSjl for the ground 
state energy is therefore adiabatic, i.e., the condition l|17|) 
corresponds to the nonadiabatic/adiabatic crossover line 
in the g/tJo, ce parameter space of the Holstein Hamilto- 
nian. Actually, the result I|18J1 is, up to the numerical 
coefficient, the same as the obtained from the symmetry 
broken side 0, |2fj, by the self-trapping of the electron 
in the continuous version of the Holstein Hamiltonian 
QJ. In this approach, the condition i|17H appears as re- 
stricting the continuous adiabatic theory to the values 
Aa -1 / 2 > 1, A < 1. For Aa" 1 / 2 ss 1 the nonadiabatic 
corrections are becoming appreciable (resulting finally in 
Eq. II12J) for small couplings). On the other hand, asso- 
ciating A -1 through A = g 2 /luqcI with the polaron width 
d by d ~ A -1 , the condition Acl keeps the continuous 
theory valid not only on the crossover line koT 1 ! 2 f=a 1 
but also for Aa -1 / 2 > 1. 

It remains thus to be shown within TPT that the in- 
frared singular (continuous), adiabatic result l|18Jl holds 
not only on the crossover line l|17JI in the 2D parameter 
space but also for Aa -1 / 2 ^> 1, as long as the continuous 
(infrared singular) approximation is valid due to A 1. 
Such a step amounts to the demonstration that TPT can 
reach beyond the nonadiabatic/adiabatic crossover itTTIl . 
associated with the symmetry breaking. Apparently this 
requires the consideration of the infinite TPT series in 
Fig. El 



g/u Q « a- 1 / 4 koT 1 ' 2 a 1 , (17) 

with a, A -c 1. In Ref.|24[ the crossover condition ifTTjl 
was derived from D'- 1 -' (at N = 1), on considering the av- 
erage number of excited phonons and on requiring that 
the latter is at most of the order of unity. Here, this 
condition follows from the electron self-energy on 
the same physical grounds, determining when the two- 
phonon processes in Fig. E|are becoming equally impor- 
tant as the single-phonon processes. 

The main contribution to £4 comes from the non-local 
(k < e + a) phonon-assisted processes. Indeed, the lo- 
cal contribution of £4^, obtained by integrating Eq. i|A2j) 
or approximately the long-wave limit of this expression 
l|A3jl over k, is apparently negligible with respect to the 
local . In other words, for e + a <C 1 the local 
contribution of the quartic non-crossing diagram and the 
non-local k m contribution of the quartic diagram are 
equally important for the values of parameters satisfying 
Eq. {JJ in determining the ground state energy and the 



IV. INFINITE SERIES 

The general diagram of Fig. can be evaluated in 
principle by using the usual diagrammatic sum rules [2lj | . 
The class of non-crossing (NC) diagrams shown in Fig. 
Q] obtained by generalizing T,f c to the order 2p, are es- 
pecially simple because they are local. The structure of 
the pth order NC diagram, proportional to g 2p , is easily 
understood by noting that the external electron phonon 
bubble differs from £2 by the fact that the single Go in 
£2 is cut in two Go's by the (p — 1) insert. Cutting Go 
in two amounts to the taking derivative of £2 with re- 
spect to a parameter in Go, analogously to the text-book 
demonstration of the Ward identities. For p = 2 this 
procedure gives immediately the result 11 -ill and iterates 
it to the order p, with /(e) given by Eq. 



The numerical coefficients A^ c are not of particular in- 
terest here, as only the scaling properties of S are con- 
sidered. 




Figure 4: The class of non-crossing diagrams corresponding 
to Eq. ltl9ll . 




q=q' 



Figure 5: Decoupling of infrared integrations over q and q' in 
Eq. (EH , 

The situation is more intricate with crossing diagrams, 
except in the particularly simple limit t « (a f=a oo), 
where the result has been obtained to all order in g. In 
order to make the discussion more transparent, k is set 
equal to zero (as appropriate for the calculation of the 
ground state energy). In addition, unlike in Eq. 119(1 or 
Eq. 1A2(I . but as in Eq. 1(16(1 . the infrared limit is taken 
from the outset on assuming 

e+(p-l)a<l, (20) 

i.e., e and a small. 

The p — 2 example of this procedure is given by Eq. 
(tl"6ll taken at k = 0. It is instructive to study first this 
example further in the limit e ~3> a (keeping e + a <C 1), 
anticipating that, for g/uiQ > oT^ the physically rele- 
vant value of e, £a, is much larger than a, i.e., that 
A ^> ojq. In this limit both contributions T,f c and 
of Eqs. l(l"5ll and iflfijl have the same leading behavior in 
e large, £4 ~ g /t 3 e 2 . This result can be derived simply 
not only for Yj^ c but for £4 too, on noting that the main 
contribution to the infrared singularity in the limit lf20)) 
with e ^> a comes from q m q 1 in Eq. 1(14(1 . Setting q = q' 
in the convolution term, the infrared integrations over q 
and q' decouple, each yielding the contribution given by 
Eq. jnj, the overall results being £^ ~ g 4 /t 3 e 2 . How- 
ever, the described approximation affects the numerical 
prefactor and the a/e corrections in £4% which are thus 



out of control in this approximation. In return, the bene- 
fit of described procedure is that it is apparently iterative 
in the sense that in the general pih order diagram it ap- 
plies, under the condition 1(2011 . to any segment of the 
type shown in Fig. This suggests that the full k = 
irreducible pih order self-energy contribution Yi 2p has the 
same form as ^ 2p c of Eq. iflfljl in the infrared limit, i.e., 



£ 2p = 2t A p (JJ -^-j a 2p (£/a) , (21) 

with the numerical coefficients A p and the function a 2p 
undetermined here. Importantly, however, for later dis- 
cussion, a 2p {£ja) tends to a finite constant for e/a large 
provided that p satisfies the condition 1(2011 . All dia- 
gramms of the sixth-order in g checked satisfy Eq. f2l")> . 

Turning to the scaling properties of £ 2p it should be 
noted that, due to its particular structure in a, the series 
ll2l)l can be rescaled from the scale 2t to the scale u)q, by 
introducing 



UJq - UJ 

i] = £ a. = 

UJQ 

to obtain 



£ 2p = c^o A p a 2p (v) (Aa^/y/r^- 1 . (22) 

This scaling of T, 2p has many satisfactory features. T, 2p is 
proportional to ujq for any p, showing explicitly that not 
only the crossing diagrams at k ~ are proportional to 
wo ("Migdal theorem") but that so are the NC ones, and 
thereby that both have to be treated on equal footing. 
In fact, the number of crossing diagrams n c is large for 
large p (n c = 1 for p = 2, n c = 8 for p = 3). Even more 
fundamentally, the proportionality to shows that the 
quantum polaron is being considered here, because ujq in 
Eq. (l22l is to be replaced everywhere by fkoo when h 
is not taken equal to 1 in Eq. QJ. Equation 122(1 also 
shows that ujq should not to be set equal to zero too 
early (or inconsistently) when the adiabatic limit is 
considered. In this respect it is also gratifying that the 
adiabatic parameter A is exhibited explicitly in Eq. 12211 . 

Equation ll22)l is also handy in the sense that, after 
using e < 1, a < 1 in the inequality 1(20(1 . to produce 
the continuous limit in Eq. 1(21(1 . it can use (the physical 
values) 77 = £/a ^> 1 to make the prefactor of (Aa -1 / 2 ) p 
small when Aa" 1 / 2 is large. This provides a strong ar- 
gument in favor of the convergence of the series 12211 . 
although the coefficients A p and a 2p (oo) are unknown. 
Here, it can be only noted that Eq. i|22|l is valid for 
arbitrary large p in the limit a < 1 provided that the 
physical relevant values of e are sufficiently small, e«l, 
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V. QUANTUM POLARON PROPERTIES FROM 

TPT 

Assuming thus that the infinite series l(22)l is conver- 
gent in the sense that it defines a function 

Z = 2taF(T 1 , Aa~ 1/2 ) , 

the physically relevant values tja of r\ are defined by Eq. 
Q). In particular, Eq. © for A can be conveniently 
rewritten as 



@. The evaluation of these from the broken symmetry 
side is quite intricate, providing the first few coefficients: 
the even coefficients Co, C2, C4, are finite and the odd ones, 
ci and C5, vanish. The coefficient c 2 is associated with 
the energy proportional to ujq and independent of A. This 
energy corresponds to the adiabatic reduction of the zero- 
point energy of the local harmonic modes which describe 
the polaron deformation. C4 corresponds to the nonadi- 
abatic kinematic corrections. The vanishing of the odd 
coefficients can be understood upon noting that when the 
leading term for A in Eq. (126)) is inserted in tja of Eq. 

(EH, 



VA - 1 = F( VA , Ka- 1 ' 2 ) , (23) 
using Eq. JHJl. This shows that va can be expressed as 



r>A - 1 = /(AqT 1/2 ) 



(24) 



Noteworthy is the fact that when £ is real, as it turns 
out to be in Eq. |(22j, the diagrammatic TPT of Eq. |23j) 
is equivalent to the Wigner TPT. The Schrodinger TPT 
leads straightforwardly to Eq. l|24*jl but through a series 
for which the convergence properties are less transparent 
than those of Eq. i'22\) . The scaling relations of the type 
(1241) are common in infrared problems [27T |. 

The function / in Eq. l(24^l has apparently two regimes 
ka- 1 / 2 1 with the crossover at ka' 1 / 2 « 1. For 
koT 1 ! 2 < 1 it can be expanded as a Taylor series in 
ka- 1 / 2 with a constant term omitted. Keeping in mind 
that A vanishes at g — this immediately gives the 
nonadiabatic result 111211 . to the lowest order in g. The 
crossover at ka' 1 / 2 w 1 is also already discussed in Eqs. 
ifTTIl and ijlSjl . as occurring towards the large, adiabatic 
Holstein polaron. The true test of the ability of TPT 
to reach beyond this crossover corresponds, as already 
mentioned, to the regime ka- 1 / 2 > 1 for A< 1. 

The only way for Eq. l(24^l to have the adiabatic solu- 
tion for koT 1 / 2 ^> 1 is that the corresponding asymptot- 
ical behavior of f{ka~ 1 / 2 ) is 



f(ka- 1/2 ) = k 2 a- 1 [c +c 1 a 1/2 k- 1 +c 2 ak- 2 +...] , (25) 



because only then a, i.e., luq ~ v M, cancels out (rather 
than goes to zero) from Eq. pi)), to the leading order in 
k~ 1 a 1 / 2 small. The final result is then 



A = 2tk 2 [c + c^^k- 1 + c 2 ak~ 



(26) 



satisfies the initial assumptions a <C 1, 
sa/o. > 1, which led to Eqs. 12112 



The result Q2 

SA < 1, VA 

Two other equivalent ways to express those conditions are 

a « 1, A « 1, ka- 1 / 2 ^ 1, or u /2t < 1, A/2i < 1, 
A/uj q > 1. 

The result Gfijlnot only reproduces the ground state 
energy (I18H [19, [2£|, but also agrees with its corrections 



A 2 a- ] 



i.e., / of Eq. 1(24^1 becomes a function of A 2 a rather 
then of ka" 1 / 2 . 

It can be objected that Eq. ll2fit follows from the as- 
sumption that TPT is convergent and adiabatic, rather 
than from the corresponding rigorous proofs. However, 
TPT is not worse in this sense than the broken symme- 
try approach. The latter also starts by assuming that 
the large adiabatic polaron exists and subsequently de- 
rives its properties and the conditions of validity. It is 
therefore all the more gratifying that those properties 
appear to be the same from the translationally invariant 
and from the broken symmetry side. 

Actually, it is possible to extend the agreement be- 
tween the two approaches from the condensation energy 
ll26)l to the polaron with d and to the polaron mass M p . 
In the continuous limit no ground state energy is asso- 
ciated with the polaron motion, i.e., the whole A is as- 
sociated with the electron localization and with the con- 
comitant displacements within the polaron width d. The 
localization energy to one site is g 2 /^o and representing 
A as A = g 2 /ujQ d, it follows that 

d = k- 1 [c' 1 +c 2 a 1 / 2 k- 1 + ...} . 

The kinematic corrections to the main Holstein result 
d ~ A -1 have not been considered before. A more rig- 
orous proof of the relation A = g 2 /ujQ d requires the 
consideration of spatial correlations, the task in which 
the diagrammatic TPT has certain technical advantage 
over the Wigner or Schrodinger TPT. 

The same type of scaling can be applied to the finite k 
states in Eq. Q and in particular to the effective polaron 
mass ©■ The derivate dT,^/duj at u — —A, associated 
with the electron spectral density, is given by the function 
F(r), ka- 1 / 2 ) of Eq. (231 . On the other hand the eval- 
uation of d£/<9£fc requires the generalization of Eq. (tl6ll . 
taken at finite fc, to all crossing diagrams. Instead of 
that it is convincing enough to note that Eq. 11 (ill makes 
it plausible to set 



1 



2t M, 



h{ka- 1/2 ) , 



p 
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in analogy with Eq. 12411 . The asymptotic behavior of the 
function h for Aa^ 1 ' 2 large must give M p linear in M ~ 
a~ 2 to the leading order in A _1 a 1 / 2 small, i.e., M v ~ 
g 8 . This also is a well-known result, derived previously 
from the broken symmetry side. On the other hand, the 
scaling properties of the electron spectral density in M, 
and then in <?, or vice-versa, remain to be determined 
because the present scaling approach only connects the 
two. 



VI. CONCLUDING REMARKS 

The central conclusion of this paper is that the prop- 
erties of the large adiabatic Holstein polaron on the ID 
lattice can be determined by the TPT summed to the 
infinite order in g. The main properties, such as the 
polaron condensation energy, its width and the effective 
mass, can all be determined from the self-energy associ- 
ated with the intermittent addition of one electron to the 
system of bosons. 

Actually, a scaling analysis rather than a precise cal- 
culation was carried out under the assumption that the 
infinite order TPT series defines a function, which has a 
meaningful asymptotic behavior in the strong-coupling 
limit: this is termed a "convergence" of TPT. The 
analysis in question was applied to the extreme limit 
ka^ 1 / 2 3> 1 and A « 1, (a « 1), which ensure the va- 
lidity of, respectively, the adiabatic and the continuous 
approximation. Although a full agreement is obtained 
with the corresponding results derived from the symme- 
try broken side it is of interest, as explained in Sec. [IJ 
to supplement the discussion with a few remarks which 
concern the applicability of TPT to other cases. 

The inequality ffity. taken in the ground state, i.e., 
with r\ = ?7a of Eq. l12.'it defines the highest order £ 2 p of 
Eqs. Q21J1 and 12211 for which A > p loq. It is obvious 
that p is the average number of phonons in the system 
excited by the added electron for a given Aa" 1 / 2 . In 
Eq. (12- >t p is taken to infinity, but it is interesting to 
consider the situation when p, i.e., Aa -1 / 2 is large but 
finite. For p < p the self-energy £ 2p is infrared singular 
and, according to Eq. |22j), £ 2p ~ (?7a) 1 ~ 3p/2 to the 
leading order in r\ large. For p < p < a -1 the self-energy 
T,2 P is still infrared singular, £ 2p ~ i r ]A)~ 3p ^ 2 (p/P)^ 1 ^ 2 ■ 
However, for p ^> a -1 (finite) £ 2p tends to the well- 
known [22] discrete, nonadiabatic t « 0, small bandwidth 
limit £ 2p ~ (1/p!) 2 . 



These observations help to elucidate the nature of the 
corrections to the large adiabatic polaron, additional to 
the non-adiabatic c 2 , C4 corrections, present already in 
the continuous limit. It appears that the Umklapp cor- 
rections to the infrared £ 2p for p < p are presumably re- 
sponsible for the build up of the PN barrier. Indeed, the 
latter are known to become important for p f» a~ x (i.e., 
for Awl), when the other candidate for Umklapp correc- 
tions, namely the infrared terms £ 2p with p < p < a" 1 , 
are squeezed out. On the other hand, the terms at 
p > a^ 1 are expected to lead to the non-adiabatic correc- 
tions to the translational dynamics of the small polaron. 
The strength of the p < p terms can be evaluated using 
a scaling procedure different from Eq. Jj^J, which shows 
that the PN potential is adiabatic, i.e., depends only on 
t and A, and is exponentially small j23| for A <C 1. 

This reasoning shows that the TPT series evolves 
smoothly from one regime to another, by changing the 
nature of the main p < p terms and the nature of its 
p ^> p tail. The present result for the large adiabatic 
Holstein polaron, together with the well-known applica- 
bility of TPT to the Lang-Firsov limit for arbitrary g, 
suggests then strongly that TPT can cover the whole 2D 
parameter space of the Holstein Hamiltonian. 

An interesting unsolved problem in which TPT can 
be useful refers to the evaluation of the polaron mass 
M p for the adiabatic polaron propagation through the 
PN barriers. Actually, the quantum tunneling through 
the PN barriers makes M p nonlinear in M in a way not 
yet clearly distinguished from that in the nonadiabatic 
regimes. 

Finally, it is clear that the problems, analogous to 
those opened and partially answered here, arise in con- 
nection to numerous quantum crossovers related to the 
symmetry breaking, and that some of the ideas devel- 
oped here may apply to those cases too. As mentioned 
in Sec. Q|this concerns in particular the problems of quan- 
tum solitons or doping the 2D Mott insulator, which are 
currently subjects of extensive theoretical investigation. 
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Appendix A: APPENDIX 

The crossing fourth-order diagram in Fig. is a convolution, 



T.^) = 9 -Y,Gk- q {^- wo)r», (Al) 



with Gfc(w) the free electron propagator and T q {ui) the leading Migdal vertex correction, the latter given by a convo- 
lution 

r 9 M = — G g /_ g (w - 2w ) G q >(uj-uj Q ) 
i' 

= ?L Z y 4 ; 

t 2 N iZ_, ( Zql _ Zq y + )( Zql - z q y-)(z q > - x+)(z q > -x-) ' 

where z q = e %q and 



7 n = 1 + na - w/2t > 1 , a;± = 71 ± yj i{ - 1 , y± = 72 ± \j 7§ - 1 , £+2:- = J/+2/- = 1 
As X-,y- < 1, by integrating over the unit circle, l/N^2 q — > §d(f)/2ir — ► $ z i =1 dz/iz2ir, one obtains 



= 72" Z <? 



3 2 / z 9 2/- 



i2 V - z q y+){z q y- ~ - z_) (x_ - z 5 y+)(a;_ - z q y-){x_ - x+) 



2 „,.„,._„.„, *"1 



£+2/+ - a; -2/ 
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i 2 (x+ - x_)(y+ - y-) (z q - 5 + )(z q -6-) 
where S± = y±/x T is real, and 6+6- = 1. Substituting T q (u>) into Eq. i|Alj) gives 

vC, ^ 5 4 -1 X+V+-X-V- I dz 



t 3 fc (x+-x-)(y+-y-) J i2v ( z - z^x+){z- z^ 1 x-){z-5+){z~6-) 
%l 



g A x\y+ - x 2 -V- 1 



t 3 (x + - X-) 2 (y + - y_) x\y+ + x 2 _y_ - 2cos(fc) 
In the continuous limit one obtains: 



x± 



l + e± \j2e + e 2 « 1 + e ± V2e 



(A2) 



y± ai l + e + a± A/2(e + a) 

E c fwl _ g 4 2Vi+V?+^ 1 

(2i) 3 e^FT^ k 2 + 2(2^+V^+^) 2 ' 1 J 

i.e., the result given by Eq. 11 fill . The contribution to the effective polaron mass in Eq. © from ££(w) in the 
numerator is obtained as 



2 (2t) 4 eVe+^ (2Vi + Ve+^) 3 t 4 A 6 g 
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